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Errors in the control of quantum systems may be classified as unitary, decoherent and incoherent. 
Unitary errors are systematic, and result in a density matrix that differs from the desired one by 
a unitary operation. Decoherent errors correspond to general completely positive superoperators, 
and can only be corrected using methods such as quantum error correction. Incoherent errors 
can also be described, on average, by completely positive superoperators, but can nevertheless be 
corrected by the application of a locally unitary operation that "refocuses" them. They are due 
to reproducible spatial or temporal variations in the system's Hamiltonian, so that information on 
the variations is encoded in the system's spatiotemporal state and can be used to correct them. 
In this paper liquid-state nuclear magnetic resonance (NMR) is used to demonstrate that such 
refocusing effects can be built directly into the control fields, where the incoherence arises from 
spatial inhomogeneities in the quantizing static magnetic field as well as the radio-frequency control 
fields themselves. Using perturbation theory, it is further shown that the eigenvalue spectrum of 
the completely positive superoperator exhibits a characteristic spread that contains information on 
the Hamiltonians' underlying distribution. 

PACS numbers: 03.67.-a; 02.30.Yy, 02.30.Zz, 82.56.Jn 
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I. INTRODUCTION 

Methods of controlling quantum^stems Q, |M ll> ^ 111 are needed to direct the course of chemical reactions @, Q , 
to determine molecular structure Qj, and to achieve quantum information processing The common goal is to 
preserve or manipulate a quantum system so that the effective evolution over a control sequence is precisely the 
desired process. The causes of unwanted dynamics include irreversible couplings to the environment (decoherence), 
static or slowly varying fluctuations in the system's Hamiltonian (incoherence), and systematic unitary errors. Here, 
we examine a common class of experimental imperfections characterized by a spatially or temporally "incoherent 
variation" in the system's Hamiltonian. In addition, we describe a method for mitigating this class of errors, and we 
report an experimental demonstration using liquid-state nuclear magnetic resonance (NMR) techniques . 

Experimental limitations present an important set of challenges to achieving precise control over quantum systems. 
Incoherent errors, which are often present in experiments, may be distinguished from decoherent because they are 
in principle refocusable. In the case of a spatially distributed ensemble interacting with a field, the field amplitude 
and direction can vary over the ensemble so that the ensemble average dynamics becomes a convex sum of spatially 
distinct unitary processes. The ensemble's dynamics will in general appear non-unitary, but as long as the correlation 
between the spins' location and the strength/direction of the field remains unchanged the dispersion in the dynamics 
can be refocused, at least in principle. 

Spatially incoherent errors have been a recurring topic of interest in the field of NMR, where they arise as inho- 
mogeneities in the static and radio-frequency (RF) fields involved. The spatially incoherent evolution caused by the 
inhomogeneities dephases the spins in the NMR ensemble, attenuating and rotating the final state away from the 
desired state. Past methods of refocusing incoherent evolution have been directed mainly towards achieving rotations 
of specific states rather than implementing a desired unitary operation on all possible states. Examples include "com- 
posite pulses" [i3j[IBIl3)[ld)llJ|j which have played an important role in creating the robust RF pulse sequences used 
for spin inversion, spin excitation, and decoupling. Adiabatic pulses have also been used to overcome the problem of 
RF inhomogeneity .15, 16J. In addition, average Hamiltonian methods have been applied to create pulse sequences 
that depend mostly on phase changes in the RF fields to control the spins, making the RF field's amplitude a less 
important experimental parameter |l7j . Most of this work has been devoted to designing robust one-spin operations, 
but more recently the work has been extended, in the context of quantum information processing, to include two-spin 
operations [l8l |. 

Here, we extend our previously described method jl^ for creating RF control gates in liquid-state NMR quantum 
information processing experiments [23, 0, l2^ l23j| . The gates use RF fields that strongly modulate a system's 
dynamics so as to refocus all the undesired evolutions and to achieve a specific desired unitary transformation with 
high precision. The method makes use of numerical optimization and complete knowledge of the system's internal 
Hamiltonian to find a modulation of the RF field such that the effective Hamiltonian over the period for which it is 
applied is precisely the desired transformation. In this paper, we extend this method so as to also use knowledge of 
the incoherent distribution of errors in an effort to design pulses that are robust over the entire distribution and thus 
yield the desired transformation with even higher fidelity. 



II. MODELING, MEASUREMENT AND ANALYSIS OF INCOHERENT PROCESSES 

In quantum information processing, the central experimental goal is to efficiently generate any desired unitary 
operation on a quantum system. Unitary operations are realized by manipulating externally-controlled parameters 
in the system's Hamiltonian. The problem of finding the appropriate external parameters, however, is an inverse 
problem, and analytical solutions are not available except for the simplest cases. As described below, we can utilize 
the more easily solved forward problem in a numerical search for a control field modulation that solves the inverse 
problem, at least in small Hilbert spaces. The density matrix evolution under an incoherent process appears to evolve 
in a larger space with a continuum of parameters, denoted by T*, describing the variations over the ensemble. In NMR, 

would be the spatial location of the individual molecules. Every measurement is an integral over T*, so assuming 
that the local density matrix p(Jr^) is uncorrelated with the local unitary operation U{~r^), the ensemble-average 
density matrix pin — J p{~f*)d~r' evolves in Hilbert space as 

Pout^ j u{-r')p,r,u\-r)dr (1) 

The variation of the operator U{~r') is tied to the variation of a field in the corresponding Hamiltonian. The evolution 
can be expressed in matrix form as a superoperator acting on Liouville space, i.e. 

\Pout) = / {u{r) u{-r)) |p,„) (2) 
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where \p) is the columnized density matrix (obtained by stacking its columns on top of each other), U is the complex 
conjugate of U and "(g)" is the Kronecker product of the matrices this form, it is clear that the input state 

\pin) is transformed by the superoperator 

s = j (u{-r)®u{-t))d-f' (3) 

to produce the final state \pout)- So although every element of the ensemble evolves unitarily, the ensemble-average 
density matrix undergoes non-unitary dynamics. 



A. Incoherent Processes in NMR Spectroscopy 

In the specific case of liquid-state NMR, the task of quantum control is to determine the RF pulse sequence that 
modulates the internal spin Hamiltonian of every molecule so as to generate a desired evolution. The homonuclear 
internal Hamiltonian for a molecule containing N spin 1 /2 nuclei is 

N N N 

k=l j>k k=l 

where — 7fc(l — a^)Bo{l^) represents the chemical shift frequency of the kth spin (7 is the gyromagnetic ratio and 
is the shielding constant), Jkj is the coupling constant between spins k and j and li denotes the i axis spin angular 
momentum operator. The chemical shifts are functions of space because the main magnetic field is never perfectly 
homogeneous througout the sample. The scalar coupling constants, however, are independent of location since they 
depend only on molecular bonding. 

The corresponding experimentally-controlled RF Hamiltonian is 

N 

k=l 

where the time-dependent functions Bfip{t) and (j){t) specify the applied RF control field, while /(~r*) reflects the RF 
field strength distribution over the sample. Here, we assume that only phase changes in time are important and drop 
the spatial dependence of (f). 

Excluding decoherence, the evolution generated by the above Hamiltonians between time and Tf is 

Ugate{~^) = Texp (^-Z^ ' [H,nt{~^) + Hext{~^,t)] d^j , (6) 

where T is the Dyson time-ordering operator. The goal is to determine a function B]^p{t) and (j)(t) that results in a 
net evolution that at every location is close to the desired transformation. In general, this inverse problem is difficult 
to solve but the forward problem of evaluating Ugate is readily solved numerically for small spin systems. 

To efficiently calculate Ugate{^), several simplifications can be made. First, the use of shim controls allows the 
static magnetic field homogeneity to be made less than 1 part in 10® over the sample volume, meaning that the 
resonance frequencies of magnetically equivalent spins differ by the same amount. As a result, for now we will drop 
the spatial dependence of Bq', the major problem is the RF field variations in H^xt- We have shown previously that 
in the fully coherent case strongly modulating RF pulses with piecewise constant RF amplitudes provide an easily 
computed modulation sequence with sufficient control over the spins . The evolution of a single pulse with constant 
(in time) RF power B^p, phase 0, frequency vpip^ and of duration r, can be solved with one frame transformation 
and a single diagonalization in that frame. This simplification allows the net propagator of a train of M such square 
pulses to be written as 

M 

Ugate = {l^RP,,n , T^m )exp [~iH"j: f {Bjip^rn , >^R.F,m Am)Tm\ (7) 

m— 1 

where U~^{i>Rp^rmTm) executes the rotating-frame transformation of the mth period and -ffe// is the effective, time- 
independent Hamiltonian in the new frame of reference 's']. Using the standard simplex search algorithm, a set of 
parameters which match Ugate to a desired transformation can often be obtained |19j . In the coherent case, we found 
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the dynamics of these pulses to be very rich, admitting many different strongly-modulating pulses for a given ideal 
unitary gate. Here we extend this solution to the case of incoherent processes in the control Hamiltonian. 

Radio- frequency power inhomogeneity was incorporated into the parameter search by tabulating a discrete his- 
togram of the RF power variations. This distribution of RF power defines a Kraus operator sum 0,|2J| instead of the 
unitary transformation of Eq. ^ , namely 

Pout = V", AkPinAl , (8) 

where Ak = ^/pkUk and pk is the fraction of the ensemble that undergoes a unitary transformation Uk- The operators 
Uk were evaluated using Eq. lO, and they differ only in the RF amplitudes involved. In NMR spectrometers the 
frequency and duration of the pulse do not vary as a function of position, and although the absolute phase does, this 
phase is unobservable here since the same RF coil is used for both transmission and reception. 

The RF inhomogeneity of our experimental setup was measured using a spin nutation experiment in which the 
transverse magnetization of the spin system was measured after on-resonance pulses of increasing duration. The 
power spectrum of nutation frequencies (Fig. ^ is a direct measure of the distribution of RF field strengths over the 
sample. 

The inclusion of RF inhomogeneity in the calculation of the fidelity, as outlined in sectioning . increases the required 
computational resources per iteration by a factor equal to the number of intervals used to define the distribution. 
Hence for the design of the pulses here only 9 intervals were used, and this distribution is superimposed on the 
experimental measurement in Fig. ^ (gray line). 



B. Metrics for the precision of coherent control 

Two useful metrics of coherent control are the state correlation and the gate fidelity . The correlation measures 
the closeness of two density matrices pideai and pout resulting from the same input state pm but evolved under different 
propagators, one of which is an ideal gate Uideai such that 

h'ldeaL dealPinUj^^^^l (9) 

while the other is the simulated non-unitary transformation given by Eq. © 

Pout = ^AkPinAl (10) 

k 

. The non-unitary transformation is trace-preserving and completely-positive. The correlation between the traceless 
parts pideai and pout of the density matrices pideai and pout is 

. - \_ trace{pideai Pout) 

'-'{Pideai, Pout) — , , (tij 

\/i™ce(p^^^^;)trace(p^„J 

which varies from —1 to 1. The correlation reports only similarities in the "direction" between two density matrices. 
To account for the loss of information due to non-unitary operations, an attenuation factor is inserted, leading to the 
attenuated correlation, which is defined by 



Ca{P ideal, Pout) — C {pideal , Pout) \ j TZZZZT^TT 

(12) 



' trace{pl^t) 
trace{pl^) 

trace{pideal Pout) 



(13) 

^/trace{pi^^Jtrace{piJ 

The attenuation term quantifies the amount of information lost due to a decrease in the magnitude of pout relative to 

Pin- _ _ _ 

The gate fidelity F is a measure of the precision of an operation 0, l25l 



F = CA{pideal,Pout), (14) 

where the overline notation Ca represents the average attenuated correlation over a complete basis of orthonormal 
Hermitian matrices pj {i.e., trace[pj pk] = 5jk). In the case of incoherent processes, it is more convenient to use the 
equivalent expression for F directly in terms of Uideai and the Kraus operators {^fe}, namely 

P = ^T.V'^^lealAk)'^ (15) 

where N is the number of qubits, which can be evaluated more rapidly than the form in Eq. (|14|l . 
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C. Exploring the Achievable Fidelities 

In this subsection, we briefly revisit the coherent case. The ideal single-spin gate fidelities of the pulses previously 
achieved values over 0.999 in numerical simulations p^ . Two important experimental parameters place restrictions 
on the achievable fidelities, the main field strength Bq, which determines how well the individual spins are resolved, 
and the maximum allowed RF power, which limits the modulation rate. To learn how changing these parameters 
affects the fidelity of a gate, we evaluted the fidelity for a 7r/2]^ alanine pulse as a function of allowed main field 
strength and RF amplitude. The results are illustrated in Fig. El 

The three curves represent the results for each of three static magnetic field strengths tested. Each point represents 
the highest fidelity found by the pulse search method out of several attempts to find the gate at that magnetic 
field and RF power. The stronger static fields enabled pulses with higher fidelities because the stronger fields cause 
the differences among the spins' frequencies to increase, allowing more lattitude for addressability and control. The 
fidelities also increased monotonically with the maximum allowed RF power. At low powers, the RF control is 
insufficient to average out the internal Hamiltonian, resulting in low fidelities, while at high RF power, the strength 
of the RF dominates the internal Hamiltonian and hence allows the dynamics to be more precisely controlled. 

It is important to emphasize that the maximum achievable fidelities of Fig. |21 represent the best gates achieved using 
the current pulse parameterization and search method and the available computational resources. These results do not 
preclude other methods and search strategies from finding higher fidelity gates. In addition, the analysis explores an 
optimistic scenario in which natural decoherence and experimental imperfections, including RF field inhomogeneity, 
are ignored. Nevertheless, expect these trends will also hold for the "self-refocusing" pulses developed in this paper. 

III. DESIGNING PULSES THAT COMPENSATE FOR RF INHOMOGENEITY 

The spatial incoherence of the RF amplitude, or RF inhomogeneity, in NMR is a consequence of the need for high 
sensitivity, which necessitates wrapping the RF coil tightly around the sample. As a result, different parts of the 
NMR sample feel unequal RF fields, causing a dephasing of the spins that attenuates the signal and introduces errors 
in the rotations of the spins. The problem can be avoided by using smaller parts of the sample (and thus less signal). 
An alternative, however, is to take advantage of the reproducibility of the field strength distribution so as to design 
gates that are insensitive to these errors. 

We used the methods of sectionmto search for both compensated and uncompensated modulation sequences for 11 
single qubit transformations. The calculations were performed for the three-spin system consisting of the ^'^C-nuclei 
in isotopically-enriched alanine. Fig. |3| summarizes the simulated fidelities for the resulting gates. 

Each fidelity point on the left of Fig. |31was calculated for a unitary gate having a single RF field strength. Each 
curve traces the gate fidelity as a function of the deviation from the ideal RF field strength. The plot shows that the 
compensated pulses are significantly less sensitive to changes in the RF modulation strength. The uncompensated 
pulses, however, have the overall highest fidelities when the RF amplitude is at its optimum value. This fact is 
confirmed in the more realistic situation where a spectrum of RF frequencies is present. The plot on the right of 
Fig. 13 graphs the pulse fidelities for the non-unitary transformations generated with RF distributions of varying 
widths. The distributions used were stretched or narrowed versions of the measured RF profile (Fig. and the 
widths have been rescaled such that a value of 1.0 corresponds to the experimentally measured RF profile. 

IV. EIGENVALUE SPECTRA OF SUPEROPERATORS 

In this section, we present a numerical study of the action of both compensated and uncompensated gates by 
looking at the eigenvalue spectrum of the superoperators of the gates. We will show how the eigenvalues of the 
superoperators can serve as a useful and convenient tool for extracting information about the imperfections in the 
implementations of unitary transformations. In particular, we shall describe features of the eigenvalue spectrum that 
allow us to characterize the distribution of incoherent processes in the superoperators, so that the reduction of these 
features in our simulated superoperators provides further evidence for the closer-to-unitary behavior of compensated 
gates as compared to uncompensated ones. 

A. Perturbation Analysis of the Eigenvalue Spectrum 

The eigenvalues of operators in SU(2^) lie on a unit circle in the complex plane. Let N be the number of spins and 
Uk denote the unitary operator determined by the RF field in the fcth frequency interval of the RF amplitude profile 
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(as shown in Fig. ^ . The eigenvalues of the superoperator obtained by taking the Kronecker product of Uk with its 
complex conjugate, tJk ® Uk, are the products of the eigenvalues of Uk with those of C/^, yielding 2^ eigenvalues that 
are equal to unity and 2^^"^ — 2^^^ pairs of complex conjugate eigenvalues. As a consequence, the weighted sum of 
the superoperators Uk (X" Uk forces the net superoperator S to have eigenvalues that are symmetric with respect to 
the real axis (which is essentially a consequence of the fact that S preserves Hermiticity) and inside the unit circle. 

Estimates of the actual eigenvalues of S will now be obtained using first-order pertubation theory. Because the 
pulses are not perfect even in the absence of RF field inhomogeneity, we may assume that the unperturbed eigenvalues 
are generically non-degenerate. The unitary operator Uk may be written in exponential form as 

Uk = e-'"" * (16) 

where Hk represents the effective Hamiltonian of the evolution over the period t for which the pulse is applied. 
Defining Hq to be the unperturbed Hamiltonian (taken to be at the peak of the profile in Fig.QJ, the eigenvalues (pj 
and eigenkets \4>j) of Hq satisfy the eigenvalue equation 

f/o|0,)=e-*^^*|0,) (17) 

where Uq = exp(— ii/oO- Note that this implies ?7o|</'j) — where \<j)j) is the complex conjugate of the ket 

\4>j)- The Hamiltonian of a particular Uk can be written as 

Hk = Hn + Kk (18) 

where Kk is the perturbation. To first order, the new eigenvalues of Hk are 

4fe = (/.j + (0j|iffe|</.j) (19) 

The corresponding eigenvalues of Uk are then 

g-i^j.fct _ g-i(<^j + {03l-Rrfc|0j>)t j-2Q^ 

Given that S = J2kPkUk <8) Uk, the spectral decomposition of 5* is 

S = T.P'^iT. ^^^'"■'^ ® E e-^^-^ 1) (21) 

km j 

and the eigenvalues of S are then given approximately by 

k 



^pi^f^~i{'l>j-4>n.)t-ii{4>j\K^\4,,)^{4>^\Kk\<p„^))t ^22) 



k 

k 



t^j„^-i{(<l>AKt,\4>,)-{4>^\Kk\4>m))t 



where in the last line we have used the fact that the diagonal elements of any Hermitian operator are real. 



B. Eigenvalue Spectra of Uncompensated Pulses 

We now use the result of the previous section, Eq. H22|l . to calculate the first-order eigenvalues of an uncompensated 
^ pulse that rotates the alanine carbon spins 1 and 2 about the x-axis. Fig. 0] compares the numerically exact 
eigenvalues with the results from the approximation. The dots are the eigenvalues corresponding to the full simulation 
of the gate under the influence of the experimental RF inhomogeneity profile. The crosses are the eigenvalues computed 
by first-order perturbation theory. To calculate Xjm, we first determined Kk using 

Kk = IX (log(C/fc) - log(C/o))A 
= Hk ~ Ho 

where log is a logarithm of the matrix obtained from the principal branch and then used Eq. (|22|l . 
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Fig.2]also shows an interesting relationship between the phase shift and the attenuation of the eigenvalues, i.e. the 
larger the phase shift, the larger the attenuation. To get more physical insight of the pattern emerging in Fig. 01 we 
now assume that can be written as — {iOk ~ ^o)K. Note that we have dropped the dependence of K with 
respect to k (i.e. the "direction" of the effective Hamiltonian is constant across the sample, which would be expected 
for a strong RF pulse). This would be in fact an exact result for a one spin system on resonance. The eigenvalue 
expression then becomes: 



Xj^ = e~ii't>3-'t>'^)t^pf^e~^^i<l>j\Kk\'Pj)-{'Pn,\Kk\<l>„^))t ^23) 

k 

— ^-i{'l>j-4>m)t'^2p^^-ii:jm(.uJk-uJo)t ^24) 
k 

thereby yielding in the continuous limit 

= e-*^'^^-'^'"'* / p{Auj)e~''''^"^^'^'dAuj (26) 



where Alo — LOk — LOo and where ipjm is an unknown real quantity. We see in this case that to first order the eigenvalue 
Xjm is just the unperturbed eigenvalue e~*'^'^J ~'^'")* times the Fourier transform of the RF distribution profile evaluated 
at ijjjm ■ This result demonstrates that the probability distribution profile of an incoherent process could be determined 
in principle from the eigenvalue structures of an experimental superoperator and a unitary-desired one. Knowing the 
form of K would indeed allow one to build the correspondence between Xjm and ipjm, and then to determine p{Auj) 
by performing an inverse Fourier transform. This result holds of course when the perturbation is in the first order 
regime and when the unperturbed eigenvalues are non-degenerate. 

To gain further understanding about the estimation of the eigenvalues of S, K^. can be expanded in the basis of 
products of Pauli spin matrices and the matrix elements contributing to the perturbation theory analysis can be 
determined. As an example, let us consider a one-spin system. In this case, can be expanded in terms of a^, (Jy 
and Uz- If we take H^t — ja^, then we see that any operator in the expansion of Kj^ that anticommutes with 
results in zero diagonal matrix elements in the eigenbasis of , because 

Aa^ + a^A^O => {±\Acr^ + a^A\±) ^ (27) 
=> ±2(±|A|±) = 

where |±) denote the eigenstates of <7x, and A is an operator in the decomposition of Kk that anticommutes with 
ax- In this simple example, only the component along Hq in Kk yields a non-zero contribution in the eigenvalue 
calculation given above (/ is not present either). In our three-spin system, if Hgt — ■ja].-, it follows that out of 
the 64 possible operators in the expansion of Kk only a]., af, a^, crfcTj, crlcFi , c^o"! and a^afaj will give a nonzero 
contribution, where i,j = x,y ot z. In general, the number of operators that yield a non-zero contribution for a 2" by 
2" unitary operator J7 is 2" — 1, i.e. the number of diagonal elements minus 1 (the identity part). This quick analysis 
reveals to some extent the slight discrepancy we have between the first order perturbation theory results and the fully 
simulated ones. Due to the nature of our physical system, the first order perturbation approach takes into account 
only a small number of elements in the decomposition of Kk- However, because it captures the general distribution 
of the eigenvalue spectrum, we can conclude that a significant amount of the perturbation is contained in these few 
operators. 



C. Symmetric Inhomogeneity Profile 

To test further our analysis and the validity of our assumptions we calculated the superoperators where a symmetric 
inhomogeneity profile is present. We recall the previous formula 

k 
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where Vjmi^^k) = {{(t>j\Kk\(f'j) ~ {4'm\Kk\(f>m)) and Aujk = uJk — (^o- Now because p{Aujk) is symmetric with respect 
to 0, one can see that the additional assumption of Vj„j(— Ao^j.) ~ —Vjjn{A(^k) (as we did in the previous subsection 
by setting Kk — AuJkK) leads to the result that Xjm — e^*^'^^"'^"-'*ylj>„ where Ajm is a real number smaller than 1. 
In other words, in the case of a symmetric inhomogeneity profile, and provided the above assumption is valid, the 
eigenvalues simply get attenuated, and are not phase shifted. As one can see in Fig. [S] for a symmetric profile, the 
eigenvalues have practically the same phase as when there is no RF inhomogeneity and are simply attenuated by the 
real factor Aj„i . This provides further evidence that the new eigenvalues are correlated with the Fourier transform of 
the inhomogeneity profile, since the Fourier transform of a symmetric real function is real. 

D. Eigenvalue Spectra of the Compensated Pulses 

As described above, pulses that are compensated for RF inhomogeneity result in an overall operation closer to 
being unitary than for uncompensated pulses. In Fig. [SJ we compare the eigenvalue spectrum of the superoperators 
corresponding to a compensated and uncompensated ^ pulse about the x axis on the alanine spins 1 and 2. 

In general, a unitary superoperator in the Zeeman basis |24| which is trace-preserving and completely positive 
must correspond to a unitary process in Hilbert space, and can be written as U i^U (Boulant et ai, in preparation). 
Since the superoperator of an incoherent process can be written as a trace-preserving Kraus operator sum, it must be 
completely positive and hence is physically reasonable. Thus the fact that the eigenvalues of the superoperators that 
were simulated for the compensated pulses basically lie on the unit circle provides further evidence for them being 
very nearly unitary. 

V. EXPERIMENTAL RESULTS 

The experimental tests were carried out on the three carbons of ^'^C-labeled alanine using a 300 MHz Bruker 
Avance spectrometer (see Ref. for the internal Hamiltonian's parameters). The experiments tested compensated 
and uncompensated versions of seven different spin-rotation pulses, which were then used to selectively evolve the two 
scalar couplings J12 and J23 while refocusing the other couplings |32|. These pulses were applied to the three input 
states 

=/j+/2 + /3, (29) 

where j = x, y, z. The input states were created using compensated or uncompensated pulses, depending on the type 
of gate being tested. For each gate tested, the average over the three input states was calculated and are shown in 
Tables III and ITU TableQlshows the results for the input states, the states obtained by selective coupling, and an average 
over all the single- pulse experiments, while Table Hll shows the specific results for each single-pulse experiment. 

The input and output density matrices were measured using state tomography and were used in Eqs. and l|12() 
to evaluate the correlation C, attenuation A, and attenuated correlation Ca {Ca — C ■ A). In all cases, the result of 
state tomography on the thermal state l\ + II ^ 1\ was used as the reference for the attenuation {pin in Eq. H12() '). 
State tomography [27j employs read-out pulses to rotate unobservable elements of the density matrix into observable 
single-spin transitions. For the three-spin alanine system used here, eight repetitions of the experiment, each with a 
different readout pulse, were used to reconstruct the density matrix |33l |. 

To determine the actual experimental gate performed, however, one requires knowledge about the effect of the 
gate on a complete set of input states so that the full superoperator can be determined, a procedure called quantum 
process tomography [23, 123, ISfll • To carry out process tomography for a single alanine gate requires state tomography 
of the input and output density matrices for a set of 64 linearly independent inputs. Thus full process tomography 
for one gate would involve 2 x 64 x 8 = 1024 separate experiments. While this is certainly possible, it is quite 
laborious. For this reason we limited ourselves to performing state tomography on just the three input states and the 
corresponding output states described above. The results do not fully characterize the experimental transformations, 
but they provide a reasonable estimate of the quality of the gates. 

The most notable difference between the results for compensated and uncompensated gates is in the severe atten- 
uation caused by RF inhomogeneity on the uncompensated gates. The average attenuation value Acomp. for the spin 
rotations is nearly four times closer to unity than the corresponding value for normal pulses (i.e. (1 — Auncomp.) / — 
Acomp.) — 3.92). In addition, Acomp. for the J23 coupling is more than 10 times closer to unity than Auncomp., although 
Acomp. for 0^12 is only slightly higher than Auncomp.- The correlations for the compensated rotation gates averaged 
0.991, slightly but consistently below the average of 0.995 for the uncompensated gates. The attenuated correlations 
for the compensated gates, however, were higher than the corresponding values for normal pulses, and the difference 
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was caused by the sharp differences in the attenuations. The results confirm that incorporating RF inhomogeneity 
into the pulse design can yield more robust gates and narrows the gap between the experimental implementation and 
the simulation results. 

VI. CONCLUSIONS 

In this report, we explored the capabilities of numerically-determined strongly- modulating pulses to achieve high 
gate fidelities in the presence of a common variety of experimental imperfections. In particular, we showed that 
robustness against incoherent errors such as RF field inhomogeneity can be obtained when knowledge of these errors 
is incorporated into the pulse design process. The basic ideas used in this paper, i.e. strong modulation, numerical 
pulse-design procedures, and the incorporation of incoherent errors in these procedures, are likely to find broad 
applicability in the development of quantum information processing devices based in many diverse physical systems. 
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TABLE I: Summary of experimental results. The metrics C, A, and Ca, refer to the correlation, attonuatiou, and 
attenuated correlation (Ca = C-A). The superscripts specify whether the pulses employed were compensated or uncompensated 
for RF inhomogeneity, while the angle brackets denote that the reported quantities are means over the three input states tested 
for each transformation. In the case of the spin-rotation values, the quantity reported is the average of all the spin-rotation 
results. 
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TABLE II: Experimental results of spin-rotation gates. The metrics C, A, and Ca, refer to the correlation, attenuation, 
and attenuated correlation {Ca = C ■ A). The superscripts specify whether the pulses employed were compensated or uncom- 
pensated for RF inhomogeneity, while the angle brackets denote that the reported quantities are means over the three input 
states tested for each transformation. The spin-rotation pulses tested were 7r/2 and tt rotations of the carbon spins denoted in 
the superscript. 
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FIG. 1: Radio-frequency inhomogeneity profile. The RF inhomogeneity in the carbon channel was measured using a 
spin nutation experiment. The resulting decaying signal was Fourier transformed to distill the various RF nutation frequencies 
present in the sample. The dotted line is the plot of the Fourier transformation, and it is the measured RF inhomogeneity 
profile. The solid gray line is the profile that was used to design pulses compensated for RF inhomogeneity, and it was extracted 
from the measured profile. 

0.99999 



0.9999 



>. 0.999 

0) 

^ 0.99 
0.9 



10^ 10^ 10^ 10* 

Maximum RF Amplitude (Hz) 

FIG. 2: Numerical exploration of achievable fidelities. The plot shows the maximum fidelities found for a 7r/2]^ alanine 
pulse as allowed RF amplitude and magnetic field strength Bo were varied. The three lines correspond to the magnetic field 
strengths explored, as denoted by the legend. 
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FIG. 3: Simulations of compensated and uncompensated pulses as a function of radio-frequency strengths and 
distribution widths. The dashed hnes correspond to compensated pulses, while the solid gray lines denote uncompensated 
pulses. The left plot shows how the compensated pulses maintain high fidelities even when the RF strength is scaled from the 
ideal value. The plot on the right simulates the same pulses as a function of the scaled width of the RF inhomogeneity profile. 
These results demonstrate the improved fidelity of the compensated pulses for all but the narrowest RF distributions. At the 
small widths, the RF profile would no longer be inhomogeneous, eliminating the need for the compensated gates. 
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FIG. 4: Eigenvalue spectrum of the simulated superoperator for a test gate and the radio-frequency inhomo- 
geneity profile shown in Fig. The dots are the exact eigenvalues of S while the crosses are the ones obtained by using 
the first order perturbation analysis. The zoom box shows some of the detail in the left-hand side of the plot. The additional 
trend line drawn in the zoomed plot emphasizes the fact that the larger the phase shift is, the larger the attenuation is. 
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FIG. 5: Eigenvalue spectrum of the simulated superoperator for the symmetric radio-frequency inhomogeneity 
profile. The dots axe the unperturbed eigenvalues and the crosses are the ones computed by using first order perturbation 
theory. The symmetry in the distribution mainly results in some attenuation of the eigenvalues with no phase shift. The zoom 
box allows more detail to be seen for the eigenvalues close to —i. 




FIG. 6: Eigenvalue spectrum of the simulated superoperators corresponding to a compensated and uncompen- 
sated pulse. The dots correspond to the uncompensated gate, while the crosses correspond to the compensated one. Note 
that the crosses basically lie on the unit circle while the dots are spread inside, confirming the closer-to-unitary behavior of the 
compensated gates. 



